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Motivated by the problem of the residual surface resistance of the superconducting radio- frequency 
(SRF) cavities, we develop a microscopic theory of the surface impedance of s-wave superconduc- 
tors with magnetic impurities. We analytically calculate the current response function and surface 
impedance for a sample with spatially uniform distribution of impurities, treating magnetic im- 
purities in the framework of the Shiba theory. The obtained general expressions hold in a wide 
range of parameter values, such as temperature, frequency, mean free path, and exchange coupling 
strength. This generahty, on the one hand, allows for direct numerical implementation of our results 
to describe experimental systems (SRF cavities, superconducting qubits) under various practically 
relevant conditions. On the other hand, explicit analytical expressions can be obtained in a num- 
ber of limiting cases, which makes possible further theoretical investigation of certain regimes. As 
a feature of key relevance to SRF cavities, we show that in the regime of "gapless superconduc- 
tivity" the surface resistance exhibits saturation at zero temperature. Our theory thus explicitly 
demonstrates that magnetic impurities, presumably contained in the oxide surface layer of the SRF 
cavities, provide a microscopic mechanism for the residual resistance. 

PACS numbers: 



I. INTRODUCTION 

Magnetic impurities in s-wave superconductors have 
been a subject of interest for a long time. Shortly after 
the development of the Bardeen-Cooper-Schrieffer (BCS) 
theory of superconductivity^, Abrikosov and Gor'kov 
(AG) demonstrated^ that magnetic impurities intro- 
duced into the sample in moderate concentrations lead 
to the suppression of superconductivity. If the magnetic 
scattering rate l/r^ > 1/r* exceeds the critical value 
1/r* « 0.88 Tco, where Tcq is the transition tempera- 
ture of a sample without magnetic impurities (we set 
h — 1 throughout the paper), the superconductivity is 
completely suppressed at all temperatures. In contrast, 
much stronger nonmagnetic disorder is required to sup- 
press superconductivity: the scattering rate 1/r ^ ep 
must be on the order of the Fermi energy ep ^ IO'^Tcq. 

According to the AG and subsequenlP^l^ theories, even 
below the critical value, l/r^ < 1/r*, the presence of 
magnetic impurities can result in the regime of "gap- 
less superconductivity" (GSC), where the superconduct- 
ing order parameter A is nonzero, yet the single-particle 
density of states (DOS) i^{e) does not vanish down to the 
Fermi level e = 0, Fig. [T] The GSC regime is predicted 
to occur quite generically, although the magnitude of the 
"subgap" DOS is parameter-dependent. Even for low 
scattering rate 1/rs ^ 1 /r* and weak exchange coupling 
J, vpJ ^ 1 {i^F is the normal state DOS at the Fermi 
level per one spin projection), optimal fluctuations in the 
impurity distribution produce "tails" in the DOS^''^ be- 
low the "hard gap" predicted by the AG theory, Fig. [ija) . 
The GSC regime becomes much more pronounced with 
increasing the exchange coupl ing and/or scattering rate. 
Fig. [ijb), as the Shiba theorj0El demonstrates. This is 
also supported by a recent numerical studjffl' of a differ- 



ent, but mathematically equivalent model. 

In the GSC regime, gapless quasiparticle excitations 
give rise to dissipation even at zero temperatur^^^. Al- 
though this dissipation mechanism (caused either by 
the natural presence of magnetic impurities or unin- 
tentional/unavoidable contamination of the sample with 
them) may be negligible for most practical applications 
of superconductors, it could play an important role in 
devices that require high quality performance. One ex- 
ample of such systems are the superconducting radio- 
frequency (SRF) cavities, widely used in particle accel- 
erators (see Ref. [T3] for a review and references therein; 
another notable system is superconducting qubits). The 
SRF cavities are characterized by exceptional quality fac- 
tors, which are, however, limited to a finite residual value 
^ lO^'^ at temperatures T <^ much smaller than the 
superconducting transition temperature Tc, where the 
contribution from thermally excited quasiparticles van- 
ishes. 

Despite the high practical relevance of the problem, 
there is no commonly accepted theoretical explanation 
of the origin of the residual Ohmic losses in SRF cavi- 
ties. Given the above properties, it was recently argued 
in Ref. [HI that they could indeed be attributed to the 
presence of magnetic impurities in the system. Although 
the bulk of Nb samples used for SRF cavities is typically 
very clean, a disordered oxide surface layer forms due 
to exposure to atmospher^Sl^ Fig. [2] Most importantly, 
magnetic moments can developP^ in the oxygen vacan- 
cies of the sub-stoichiometric Nb205 layer of thickness 
~ 5 — lOnm. Since the typical penetration depth of the 
electro-magnetic field is ^ 45nm, this means that super- 
conductivity could be partially suppressed in a consider- 
able fraction of the operating region. Scanning tunnel- 
ing spectroscopy measurements of the SRF cavity sam- 
ples performed in Ref. Uni revealed the surface tunneling 
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FIG. 1: (Color online) Regimes of "gapped" [(a), (c), (e), weaker exchange coupling ufJ and lower scattering rate 1/ts\ 
parameters of the Shiba theory used: l/ijs/S.) = 0.04, 7 = 0.95] and "gapless" [(b), (d), (f), stronger exchange coupling and/or 
higher scattering rate; parameters used: 1/{tsI\) = 0.17, 7 = 0] superconductivity, (a) and (b) The single-particle density of 
states (DOS) v{t), obtained from the Shiba theory Eqs. ( |4.8[ )-( [47l2| , and ( |6.1[ ). The inset in (a) schematically shows (in red) 
the "tail" of the DOS produced by optimal fluctuations of impurity distribution - exponentially small nonperturbative e ffect , 
not captured by the AG-Shiba theory and not considered in the present paper, (c) and (d) The function Q2(e, fc) [Eq. (5.3l] 
describing the dissipative contribution to the current response from a given energy e, see Sees, [y] and |VI| for details. The plots 
are presented for fe = 0; the inset in (d) shows the full range of Q2(e, fc). The function Q2{e,k) is nonzero if and only if i^(e) 
is nonzero, (e) and (f) The temperature dependence of the surface resistance R{uj), obtained from the main Eqs. ( |3.1[ ), ( |3.3[ l, 
(4.11, (4.151, (4.161, and (4.171 of the paper by numerically calculating the integrals over e and k. In the gapped regime (e), 
. ?([j) oc exp(— A*/r) is exponential at lower temperatures and vanishes at T = 0. In the gapless regime (f), for moderate 
"subgap" DOS, the surface resistance -R(a;) is exponential at lower but finite temperatures and saturates to a nonzero value at 
T = 0. The latter case reproduces the commonly observed experimental behavior^. 



DOS with appreciable "subgap" contribution, consider- 
ably greater than one would expect from a high-purity 
Nb material. Combined with good fits to the Shiba the- 
ory, these data suggested magnetic impurities in the ox- 
ide surface layer as an important contributing factor to 
the dissipation in SRF cavities. 

To support this idea, in the present work, we develop 
a microscopic theory of the surface impedance of s-wave 
superconductors with magnetic impurities. According to 
the surface chemistry of the air-exposed Nb samples^. 
Fig. [2] the real SRF cavity material is most appropriately 



described by a model of disordered surface layer that con- 
tains both magnetic and nonmagnetic impurities, while 
the rest of the sample is weakly disordered or pure. In 
principle, such model can be studied in the framework of 
the quasiclassical approach to s uperconductivity based 
on the Eilenberger equatiorF^'^. However, this involves 
solving a self-consistency problem for a system of differ- 
ential equations, which, for realistic parameter values, is 
challenging even using numerical methods. 

Instead, here we consider a simpler model of a super- 
conducting sample with uniform in space distribution of 
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FIG. 2: (Color online) Typical structure of the surface layer 
of the air-exposed Nb samples used for SRF cavities. Lo- 
calized magnetic moments (shown as green arrows) can form 
(Rcf. 15) in oxygen vacancies of the Nb205 layer of thickness 
^S-lOnm. Typical penetration depth of the radio-frequency 
electromagnetic field is ~45nm. 



magnetic and nonmagnetic impurities, Fig. [3] The main 
practical advantage is that for this model we are able to 
analytically obtain the general expressions for the current 
response function and surface impedance. The expres- 
sions are valid, within the approximations of the theory, 
in a wide range of parameter values and, in the general 
case, only the resulting integrals need to be calculated 
numerically. This generality allows for the application 
of our theory to the description of experimental systems, 
such as SRF cavities and superconducting qubits, in var- 
ious practically relevant regimes. On the other hand, if 
necessary, explicit analytical expressions can be obtained 
in numerous limiting cases. 

The current response function and surface impedance 
of superconductors without magnetic impurities are pro- 
vided by the Mattis-Bardeen^^- and Abrikosov-Gor'kov- 
Khalatniko\^2SI theories (see also Ref. I21)22p , in the pres- 
ence and absence of nonmagnetic disorder, respectively. 
The surface impedance of superconductors with magnetic 
impurities was previously studied in Refs. I23I24I in the 
limit lypj <^ 1 of weak exchange coupling of the AG 
theory. In the present work, we consider the case of 
arbitrary exchange coupling strength vpj, treating the 
interactions of conduction electrons with magnetic im- 



purities within the framework of the Shiba theory. This 
sufficiently widens the range of impurity concentrations, 
where the GSC regime (of particular interest to us) with 
appreciable DOS at the Fermi level occurs. 

As the feature of key relevance to SRF cavities, we 
demonstrate that the presence of magnetic impurities 
does lead to the saturation of the surface resistance at 
zero temperature in the GSC regime. 

Our theory employs the linear response formalism and 
is therefore valid as long as the superconducting state 
is not appreciably suppressed by the magnetic field H. 
For type-II superconductors in thermal equilibrium, the 
upper bound for this is set by the first critical field Hd, 
above which the system becomes unstable towards cre- 
ation of vortices. Real SRF cavities, however, are known 
to operate in a metastable vortex-free state that per- 
sists up to a higher "superheating" field iJgh > HcP^^^- 
Thus, our theory should be applicable in the range H < 

At higher fields H ^ Hsh one could, in fact, expect 
a cooperative effect of the two dissipation mechanisms: 
magnetic disorder could create "hot spot" regions^^ of 
locally suppressed superconductivity at the surface and 
thus trigger proliferation of vortices. Such regime de- 
serves a separate study. 

The rest of the paper is organized as follows. In Sec. [Ill 
the studied system is presented and the main approxima- 



tions are formulated. In Sec. Ill the surface impedance 
and current response function are introduced. In Sec. |IV[ 
the current response function is calculated. In Sec.|Vj the 
low- frequency expansion is performed. In Sec. |VI[ the key 
result pertaining to the presence of magnetic impurities 
- finite residual surface resistance in the GSC regime - 
is demonstrated. Concluding remarks are presented in 
Sec. En 



II. MODEL 

We assume the superconducting sample occupies the 
half-space z > and contains both nonmagnetic and 
magnetic impurities, which are distributed uniformly 
in space with densities n and tIj, resp ectiv ely. Fig. [3] 
Within the framework of the BCS theorjMiSfll, the Hamil- 
tonian of the system can be written as 



H 



d^r 



z>0 



r 



Here, -00- = "0,7 (r) is the electron field operator, a, a' =t, i 
are the spin indices, cr = (cr^;, Cj,, (Tz) is the vector of Pauli 
matrices, and summation over repeated spin indices is 
implied; E(p) is the electron spectrum, which we assume 



isotropic in momentum p, p = |p|, p = — iV; A = A{t, r) 
is the vector potential of the electromagnetic field pene- 
trating the sample; A is the superconducting order pa- 
rameter, which has to be found self-consistently in the 
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FIG. 3: (Color online) Studied system: a half-infinite {z > 0) 
superconducting sample with magnetic (green) and nonmag- 
netic (red) impurities uniformly distributed over its volume. 
Relation between the wave-vector k, electric E and magnetic 
H fields, and the vector potential A is shown. 



presence of magnetic impurities. 

Next, Ta and r;, are random positions of nonmagnetic 
and magnetic impurities, respectively. We use the con- 
ventional disorder averaging technique's) ( "noncrossing" 
approximation) and, to keep calculations simpler, assume 
contact interaction potential of impurities ("point disor- 
der"). We (i) treat magnetic impurities as classical spins 
described by the unit vectors and assume them unpo- 
larized, (ii) consider arbitrary exchange coupling strength 
vpJj summing the full perturbation series for a single im- 
purity. These are the approximations of the Shiba the- 
orjEl. 

Note that the exponentially small subgap contribution 
("tail") to the DOS arising from the optimal fluctuations 
of magnetic disordei^^tiSl jg not captured within the non- 
crossing approximation. This nonperturbative effect is 
dominant only in the limit of weak exchange coupling 
i^pJ ^ 1 and small impurity concentration, such that 
^/{tsTco) ^ 1- In this work, we concentrate on more sig- 
nificant contributions to the DOS that arise at larger im- 
purity concentration and/or stronger exchange coupling. 



III. SURFACE IMPEDANCE 

The complex surface impedance 

Z{uj) = Riuj) +iX{uj) (3.1) 

(see, e.g., Refs. I22l3l|) relates the electric E(2;)e~''^* and 
magnetic H(z)e~"^* fields at frequency oj at the interface 
z = between the vacuum and sample as 



E(0) = -Z(c.)[H(0) X no], 



(3.2) 



where ng = (0, 0, 1) is the unit vector normal to the sur- 
face pointing into the sample. Fig. [3 The real part R{uj) 
of the impedance ( |3.1[ ) determines the energy flux (aver- 
aged over the oscillation period) 



of the electro-magnetic fleld per unit area from the vac- 
uum into the sample and is referred to as the "surface 
resistance" . 

The generaP^l expression for the surface impedance of 
a uniformly disordered system reads 



Z{uj) = -1 5 , A(w) 



+ 00 



dk 



■AT:Q{oj,k)/c^' 
. (^-^^ 

Here, Q(w, k) is the linear current response function of 
an infinite sample in Fourier representation, dependent 
on the frequency u and the absolute value k — |k| of the 
wave- vector k. In Eq. (3.3), we also introduce the com- 



plex penetration depth A(cj): its real part ReA(a;) is the 
actual penetration depth that determines the decay scale 
of the electromagnetic field into the bulk. The Ohmic 
dissipation is determined by the imaginary part (32 (w, k) 
of the current response function 



Q{(jj, k) = Qi(w, k) - iQ2(w, fc)- 



(3.4) 



According to Eq. (3.3), the surface resistance R{uj) is 
finite, only if Q2(w, k) is nonzero. 

The response function Q[u)^k) defines the relation 



JHk) 



-Q(w,fc)A(w,k) 



(3.5) 



in the Fourier representation between the electric cur- 
rent j(w, k) and the electro-magnetic field, described by 
the vector potential A((u;,k). The trivial tensor struc- 
ture of Eq. (3.5 1 holds for cubic crystal symmetry and, 



in particular, for the isotropic electron spectrum i?(|p|) 
assumed here. It is convenient to work in the gauge of 
absent scalar potential 'p{uj, k) = 0. Additionally, the 
vector potential A(cj,k) may be assumed to satisfy the 
constraint 



A(w,k)k = 0, 



(3.6) 



which significantly simplifies the calculations. This is 
equivalent to the local electroneutrality condition, which 
is a very good approximation for superconductors. The 
geometric relation between the wave-vector k, vector po- 
tential A(aj,k), and electric E(cL;,k) = i^A(tj,k) and 
magnetic H(ci;, k) = [ik x A(a;, k)] fields is shown in 
Fig.H 

In the next section, we calculate the current response 
function Q{Lu,k), which fully determines the surface 
impedance (3.3 1. 



IV. CURRENT RESPONSE FUNCTION 



According to the general Kubo formalisrrPSllil^ the 
expression for the linear current response function 



Eq. (3.5)] can be written down as 



5 



Qiuj,k) = Qo- 



-3i 



+ 00 



de 



tanh 



2T 



tanh (jj)^^(e,.;, k) + tanh ^ Oj)««(e, k) - tanh |±;0j)^^(6, ^, A:) 

(4.1) 



"a / ^[G"(6+,p+)G^(e_,p_)+F"(6+,p+)F^(e_,p_)]^ , a, 6 = i?, A 



(4.2) 



In Eq. (4.1), we introduced 



(4.3) 



where vp = idE/dp)p=pp is the Fermi velocity, is the 
Fermi momentum, E(pp) = ep, and = p'p/i'^T^'^vp). 
The quantity Qo = <5(0, 0)|t=o, clean is the value of the 
current response function for a clean system at w = 0, 
k = 0, and T = 0; it can be related to the formal Lon- 
don penetration depth Alq (which can be introduced as a 
characterization parameter, regardless whether the Lon- 
don limit actually applies) as 



1/Aio - 47rQo/c^ 



(4.4) 



The quantities Qq and A^o describe the corresponding 
clean system and are determined only by the band struc- 
ture parameters. 



In Eq. (4.2) 



A 



g2 _ ^2 _ ^2 

G^{e,p) = [G^(e,p)]*, F^{e,p) = [F«(e,p)]* 



g2 _ ^2 „ ^2 ' 

(4.5) 



are the retarded (R) and advanced (A) "normal" (G) and 
"anomalous" (F) Green's functions, averaged over disor- 
der, where the conventionally introduceci^ functions e 
and A are de fined below. For point disorder and due to 
the property (3.6 1, the "ladder" contributiorP^ vanishes. 



and the current-current correlation functions (4.2) are 



determined by the products of disorder-averaged Green's 
functions. 



Further, in Eqs. (4.1), (|4.2[), and (4.5 1, e is the energy 



relative to the Fermi level ep, and ej- = e ± uj/2, p± 
p ± k/2. We split the integration over momentum 



d^p 

(27r)3 



dn 

|n| = l 47r 



in a standard way into the integration over its absolute 
value p = IpI, expressed in terms oi ^ = E{p) — ep, and 
averaging over its direction, expressed in terms of the 
unit vector n = p /p. 

Equation (4.2) defines the correlation functions of the 
current components perpendicular to k, see Eq. ( |3.6[ ), 
and so, n„ are the components of n = {nx,ny,nz) per- 
pendicular to k: if k = (0,0, A:), as in Fig. [s] then 



a = x,y. Since the spe ctru m is isotropic, the integral 
/(nk) — / ^ . . . in Eq. (4.2 1 depends just on nk = n^k 
and angular averaging is reduced to calculating the inte- 
gral 



dn^{l - nl)I{n^k). 



(4.6) 



Note that for the integration order as in Eqs. (4.1 1 and 



(4.2) - first over ^ and then over e - the contribution 



to Q{uj,k) arising from the dependence of the current 
operator on the vector potential is already compensated 
for. 

It is convenient to i ntrod uce the (retarded) quasiclas- 
sical Green's function^^l^ 



(4.7) 



{5,/}(e) = ^ J daG",F"}{e,p), 



which, according to Eq. (4.5), equal 



Ve^ - A2 Vv^ 



fie) 



- A2 Vv^ - 1 ' 
Within the Shiba theory^, the function 
. X _ e 5(e) 

^ = "^^^ = A-7(^ 

satisfies the equation 



vA = e- 



1 Vl - u2 
Tg 7^ — ^2 



Here, 



7 



1 - (Triypjy 
1 + {iTiypjy 



(4.8) 



(4.9) 



(4.10) 



(4.11) 



(4.12) 



is the parameter of the Shiba theory characterizing ex- 
change coupling strength and is the scattering time on 
magnetic impurities, 

'^^ -(1 - 7^) = 27r^'Fns 



T, 2'Kl'p 



[1 + {niypjy 



6 



In the weak coupling limit lypJ <C 1 of the AG theorj-™ In Eq. (4.131, t is the scattering time on nonmagnetic 
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We also introduce the function h{e) as 



Ve^ - A2 = hie) + ^, 
2t 



which is related to g(e) and /(e) as 
1 / e A 



He) = 



2 V.9(e) fie) J ' 



(4.13) 



(4.14) 



impurities, 



= 2'KVpnu 



Solving Eq. (4.11 ) for v (in the general case - numeri- 
cally), o ne ob tains g{e), /(e), and hie). Integration over 
^ in Eq. (4.2) is straightforward and we obtain 



(jj)«^(e,a.,fc) = ^[g(e+)g*(e_)+/(e+)/*(e_)+l](jj)«^(e,..,fc), (jj)«^(e, fc) = ( ;,(,^) _ fe. (/"j"^ , ^^^i, 

(4.15) 

(jj)««(e,^,fc) = i[l-g(e+)g(e_)-/(e+)/(e_)](J7)««(e,..,fc), (j7)^'«(e, fc) = ^ ^ , , "f"., , / , 

2 \ /i(e+) + ft(e_) + i/r - i;Fnk/ ^ 

(4.16) 



r 



and {jj)^^ie,uj,k) = [(jj)^«(e, A:)]*, (jj-)^^(e^fc) 
\{ij)^ ^ie, LP, fc)]* . Angular averaging in Eq. (4.151 and 



(4.16) can also be performed explicitly according to 
Eq. (4.61), 



(JJ) 



RR,RA 



ie,u,k) = 



«0 + -L 



(4.17) 



where 



to = — r X 



/i(e+) - ft,*(e_) + i/r, for 
7i(e+) + /i(e_) + i/T, for 



Equatio ns ([4l|), (|4.15| ), ( |4.16| , and ( |4.17p , combined 
with Eqs. (4.8)-(4.14) of the Shiba theory, provide the an- 



swer for the current response function Q(w, k) and consti- 
tute the main result of our work. Within the approxima- 
tions of the theory, these equations are valid at arbitrary 
values of frequency oj, temperature T, and six micro- 
scopic paramete rs cha racterizing the system. Three stan- 
dard parameterJlSHni describe the clean system: (i) the 
superconducting transition temperature T^o or, equiva- 
lently, the superconducting order parameter Ap at T = 
for a system without magnetic impurities; (ii) the current 
response Qo [Eq. (4.3)] or, equivalently, the formally in- 
troduced London penetration depth Alq [Eq. ( |4.4[ )] for a 
clean system at T = 0; and (iii) the Fermi velocity vp. 
The other three parameters describe disorder: (i) the 
nonmagnetic scattering time r; (ii) the magnetic scatter- 
ing time Ts] and (iii) the exchange coupling strength vpj 
or, equivalently, the Shiba parameter 7 [Eq. (4.12)]. 

In the general case, Eqs. (4.1), (4.15), (4.16), and 
(4.17) provide the most explicit analytical form of the 
current response function (3(w, k) possible. The func- 



tion is given by the integral over energy e in Eq. (4.1), 



value of momentum k is explicit in Eqs. (4.17), while 
the dependence on e is obtained from the well-known 



Shiba equation (4.11), the solution to which determines 



where the dependence of the integrand on the absolute 



the functions g(e), /(e), and hie). 

For arbitrary values of parameters, the solution to the 
Shiba equation and the integrations over e for the current 
response Q(w, fc) and over k for the surface impedance 
Ziuj) [Eq. ( |3.3[ )] need to be carried out numerically. The 
generality of the obtained results, however, should make 
them applicable to a variety of realistic experimental 
regimes. 

On the other hand, in a number of limiting cases, the 
general formulas can be further simplified and in many 
cases explicit analytical expressions for the current re- 
sponse function and surface impedance can be obtained. 
Analysis of such limiting cases is straightforward and we 
do not present it here. 

In the weak coupling limit vpJ ^ 1 of the AG the- 
ory^ (7 = 1), the results of Refs. l23l2H are (presumably) 
recovered, and in the complete absence of magnetic im- 
purities [l/r, = 0, hie) = Ve^ - A2, g{e) = e/He), 
fie) ~ A//i(e)] the Mattis-Bardeen theorylSl jg repro- 
duced. 

In the next two sections, we consider the low-frequency 
limit io <^ A, most relevant for practical applications 
to SRF cavities, and concentrate on the key property 
pertaining to the presence of magnetic impurities - finite 
residual surface resistance in the GSC regime. 



V. LOW FREQUENCY EXPANSION 

The typical operating frequencies of the SRF cavities 
are w ~ c/L ~ 10~2meV <C A '-^ ImeV, where L ^ 10cm 
is the typical size of the cavity. Thus for practical applica- 
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tions to SRF cavities, one may perform the low-frequency 
expansion in w <C A. Separating the real ( "nondissipa- 
tive") and imaginary ( "dissipative" ) parts [Eq. (3.4|] of 



Q{uj, k) [Eq. (4.1 1], in the leading order in uj for each, we 
obtain, 



3i 



+00 



detanh— {[/(e)]20j)««(,,0,fc) 



[r(^r(.7jy^(e,o,fc)}. 



(5.1) 



+ 00 



de - 



dno(e) 



de 



(52(6, fc). 



(5.2) 



[r(£)]'(jjr^(^,0,fc) + [1 + |5(e)P + \f{e)\']{jj)^^{e,0,k)} 



(5.3) 



Here, no(e) — l/[exp(e/r) + 1] is the Fermi distribution 
function. 

The real part Qi{k) [Eq. (5.1)] is finite at w = and 



determines the penetration depth A(a; = 0) [Eq. (3.3)] of 
the quasistatic magnetic field (Meissner effect). On the 



other han d, the imaginary part Q2(u},k) (x w [Eqs. (5.2) 
and (5.3)], which determines the dissipation, is nonzero 
only at finite frequency and is linear in ilP^at I/ts 3> w. 
Since Q2{'^, k) is smaller than Qi{k) at least in w/A, one 
may also expand Eq. (3.3) in Q2{uj,k) to obtain 



32 71UJ 



r-\-oo 



dk- 



Q2(w,fc) 



[fc2 + 47rQi(A:)/c2 



(5.4) 



Thus the surface resistance R{u!) cx is quadratic in fre- 
quency, which is the most common dependence observed 
experimentally in SRF cavitieJl^. 



VI. RESIDUAL SURFACE RESISTANCE 

We now turn to th e ke y finding of our work. The func- 
tion Q2{e,k) [Eqs. (5.3|] describ es th e contribution to 
the dissipative part Q2\^, k) [Eq. ( |5.2[ )] of the current re- 
sponse function from quasiparticles with a given energy 
e, while the derivative — dno(e)/de constrains their distri- 
bution to t he ra nge |e| ^ T around the Fermi level. As is 
well knowiP^Ml^ y-eal part gi{e) of the normal Green's 
function (7(e) = 51(e) — ig2{^) [Eq. (4.8 1] determines the 
DOS 



v{e) = VFgi{e)- 



Inspecting Eqs. (5.2) and (5.3), we notice that 



'h{t,k)={)^v{e)^0. 



(6.1) 



(6.2) 



Indeed, if 51(e) = 0, i.e., g(e) = ~iff2(e) is imaginary. 



then, according to Eqs. (|4.8[), (4.9), and (4.14), so are 



/(e) = — i/2(e) = — i-\/l + .gf (e) and h{t). In this case, 
(j7)''^(e,0,fc) - (lj)««(e,0,fc) = (jI)-^-"(e,0,fc) and the 
function (|5.3|) does vanish. 



02(e,fc) = ^(jj) 



RA 



(e,0,fc)[-2/|(e)+l+g2(e)+/|(e)] = 0. 



We do not present a more cumbersome rigorous proof of 
the converse here. Instead, the property (6.2) is clearly 
illustrated in Figs, [ija), (b), (c), (d). 

Thus, as one would intuitively expect, only the energies 
e at which the DOS v{e) is nonzero contribute to dissipa- 
tion. At T = the envelope function — dno(e)/de — )■ (5(e) 
becomes a delta-function and only the excitations at the 
Fermi level e = contribute. 



j(w, fc)|T=o = Qo'^'92(e = 0, k). 



(6.3) 



According to (6.2), as the central result, we obtain that 



the system exhibits finite surface resistance (5.4) at T = 



if and only if the DOS at the Fermi level is nonvanishing, 
i.e. the system is in the GSC regime. 

R{uj)\t=o > ^ zy(e = 0) > 0. 

The result is illustrated in Figs.jl] If the spectrum has 
a gap A (DOS "tails'^ISHini are not captured by the AG- 
Shiba theory). Figs. [T]^a),(c),(e), the surface resistance 
obeys an activation law R{u}) oc e~^/'^ at temperatures 
T ^ A, eventually vanishing at T = 0. 

On the other hand, in the gapless regime. 
Figs. [T|b),(d),(f), the surface resistance i?(aj) satu- 
rates to a finite value at T = 0. We note that in 
experiments it is quite typicaP^ for i?(w) to exhibit 
both the saturation at T = and an activation behavior 
i?(w) (X exp(— A*/T) at finite but low temperatures 
T < A* , with A* close to the value of the superconduct- 
ing order parameter A of a clean sample. Such behavior 
can be realized, if the finite DOS v{e = 0) <C j^f at the 
Fermi level is much smaller than the DOS v{e > A) ^ i/p 
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above the "nominal" gap. In the framework of the Shiba 
theory, this is possible in the limit of low magnetic 
scattering rate I/t^ ^ A and stronger exchange cou- 
pling vpJ ^ 1, the case shown in Figs. [l];b),(d),(f). 
Thus, our microscopic model can reproduce the typical 
experimental temperature dependence of the surface 
resistance. 



VII. CONCLUSION 

In conclusion, we developed a microscopic analytical 
theory of the surface impedance of s-wave superconduc- 
tors with magnetic impurities. The theory can poten- 
tially be applied to a variety of superconducting systems 
and is of direct relevance to the problem of the residual 
surface resistance of SRF cavities. We explicitly demon- 
strated that, in the regime of gapless superconductivity, 
the system exhibits saturation of the surface resistance at 
zero temperature - a routinely observed, but largely un- 
explained experimental feature. This substantiates the 
recent conjectural^ that magnetic impurities, formed at 
the surface of the oxide surface layer, could be the dom- 
inant dissipation mechanism limiting the performance of 



the SRF cavities. Our theory is valid in the wide range of 
parameter values and can be used for direct comparison 
with experimental data, as will be presented elsewher^^H 
Note added. After a prep rinlP^ of the present work be- 
came available, a papep^S came out, in which the same 
problem was studied in the diffusive limit (rA ^1) and 
for weak exchange coupling {vfJ ^ 1), with an empha- 
sis on the quasiparticle contribution from the vicinity of 
the spectral gap. The approach of Ref. [3HI based on the 
Keldysh formalism and Usadel equation is equivalent to 
ours in that limit and most of the analytical results of 
Ref. [38] can be obtained as asymptotic expansions of our 
general formulas. 
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